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Abstract. We introduce a force-directed algorithm, called Sync-and- 
Burst, which falls into the category of classical force-directed graph draw¬ 
ing algorithms. A distinct feature in Sync-and-Burst is the use of sim¬ 
plified forces of attraction and repulsion whose magnitude does not de¬ 
pend on the distance between vertices. Instead, magnitudes are uniform 
throughout the graph at each iteration and monotonically increase as 
the number of iterations grows. The Sync-and-Burst layouts are always 
circular in shape with relatively even distribution of vertices throughout 
the drawing area. We demonstrate that aesthetically pleasing layouts are 
achieved in 0(n) iterations. 


1 Introduction 

Force-directed graph drawing algorithms have been established as probably the 
most successful approach towards drawing generic undirected graphs [5]. While 
the classical force-directed algorithms, such as the spring embedder of Fades |4j, 
the algorithm of Fruchterman and Reingold [5] , the barycentric approach mm 
and the algorithm of Kamada and Kawai [7] perform well on relatively small 
graphs, it has been observed that they do not scale up well. The larger the input 
graph is, the more likely a force-directed algorithm is to get trapped into a local 
minimum and produce a poor result [8]. There has also been recent work on 
improving force-directed layouts by introducing a third force in addition to the 
forces of attraction and repulsion between vertices [T]. 

The algorithm introduced in this paper is a force-directed algorithm which 
follows the general scheme of the spring embedder algorithm. However, we use 
simplified forces of attraction and repulsion which do not depend on the distance 
between vertices. Instead, both the magnitude of the attraction force and the 
magnitude of the repulsion force are uniform throughout the graph at each given 
iteration, and they monotonically increase (attraction slower than repulsion) as 
the number of iterations grows. Unlike in other force-directed algorithms, we do 
not aim at minimising the total energy of the system. Instead, we experimentally 
demonstrate that aesthetically pleasing layouts are achieved in 0{n) iterations. 
This approach results in straight-line layouts with a circular shape and generally 
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with a more even distribution of vertices throughout the drawing area when com¬ 
pared to the results of the algorithm of Fruchterman and Reingold. Furthermore, 
we show some evidence that our algorithm scales up better than the algorithm of 
Fruchterman and Reingold and is able to find highly symmetrical layouts of some 
classes of graphs which do not have symmetrical Fruchterman-Reingold layouts. 
These are graph whose symmetrical layouts would require relatively long edges, 
e.g. some Hamiltonian graphs including queen graphs, the Wagner graph and 
the Heawood graph la¬ 
in section we introduce the basic terminology necessary for presenting our 
algorithm. The algorithm is described in detail in section Section presents a 
comparison between our layouts and Fruchterman-Reingold layouts of the Rome 
Graphs and a few other selected graphs. Finally, in sectionj^we draw conclusions 
from this work. 

2 Preliminaries 

In this work we consider undirected graphs. An undirected graph G(V, E) consists 
of a set of n > 0 vertices V = {vi : 1 < i < n} and a set of m > 0 edges 
E = {ci : I < z < m}, such that each edge is an unordered pair of vertices. Two 
vertices Vi,Vj G V are adjacent if {vi,Vj} G E. Two edges are adjacent if they 
share an vertex. The adjacency matrix A = {a)ij of G is an n x n (0, l)-matrix 
such that Gij = I iff {vi,Vj} € E. 

For simplicity, we refer to undirected graphs simply as graphs in the remain¬ 
der of this paper. By a layout of a graph we understand a placement, i.e. x- and 
y-coordinates, of all vertices within a given 2D drawing area. Assuming forces 
of attraction and repulsion are defined between pairs of vertices and considering 
the graph as a mechanical system, a classical force-directed algorithm starts with 
a random initial layout and then iteratively finds a layout with a minimal total 
energy, i.e. a layout that corresponds to a mechanical equilibrium of the system. 
We denote a particular iteration of the algorithm by t with t G [1, -boo]. 

Generally speaking, a force is a vector with a magnitude and direction. We 
denote the total force of attraction between vertex Vi and its adjacent vertices at 
iteration t of a force-directed algorithm by fa(vi,t)- Similarly, fr{vi,t) denotes 
the total force of repulsion between Vi and other vertices at iteration t. We use 
fa(,Vi,t)-x, faivi,t)-y^ frivijt)-^, fr{vi,t)-y to denote the X- and y-components of 
faivi,t) and fr(vi,t) in 2D Gartesian space, respectively. Similarly, Vi(t).x and 
Vi(t).y represent the x- and y-coordinates of vertex Vi at iteration t, respectively. 

3 The Sync-and-Burst Algorithm 

The algorithm that we propose consists of two phases. The first sync phase 
brings highly interconnected vertices close to each other, while the second burst 
phase spaces vertices out. Thus, we call our algorithm Sync-and-Burst. It is 
roughly based on the algorithm of Fruchterman and Reingold [S] , but with some 
very significant differences. In the first two parts of this section we introduce 


the forces of attraction and repulsion used in Sync-and-Burst and then we give 
details about the method we use for determining the x- and y-coordinates at 
each iteration. The complete algorithm is outlined in section |3.4[ 


3.1 Forces and Total Magnitude 

Similar to Fruchterman-Reingold, we consider forces of attraction between adja¬ 
cent vertices and forces of repulsion between any pair of vertices. However, at any 
particular iteration t of Sync-and-Burst, the repulsion between any two vertices 
has the same uniform magnitude M{t) > 0 which is monotonically increasing 
with the growth of t. Similarly, at any particular iteration t, the attraction force 
between any pair of adjacent vertices has the same magnitude which monoton¬ 
ically increases as the number of iterations grows. That is, in comparison to 
Fruchterman-Reingold we take a simplified approach where the magnitudes of 
the attraction and the repulsion forces between vertices do not depend on the 
distance between them. 

We want attraction to be the predominant force in the system in the first 
iterations, making adjacent vertices sync their positions with each other in the 
layout. Then, as the number of iterations grows, we want the repulsion force to 
catch up and make vertices go away from each other at an increasing rate. This 
can be achieved by having the magnitude of the attraction force between a pair of 
vertices grow slower than the magnitude of the repulsion force. In this work, we 
set the magnitude of the attraction force to ® (however, the parameter 

0.9 can potentially be tuned to a different value through experimentation). 

With this choice of magnitudes, we consider the total magnitude M.a{t) of 
attraction and the total magnitude of repulsion A4r{t) in the system (not the 
magnitude of the vector sum of forces) as an estimate of how much attraction 
and how much repulsion there is in the system at a particular iteration t. These 


total magnitudes and their difference f{t) = Ma{t) — Mr (t) 
Equations Q-([^. We call f{t) the total magnitude function. 

are expressed in 

Ma{t) = 2M{tf-^m'^ 

(1) 

Mrit) = M{t)n{n — 1) 

(2) 

f{t) = 2M(t)°®m^ — M{t)n{n — 1) 

(3) 


A typical behaviour of f{t) for a monotonically increasing M{t) is plotted 
in Fig. The part of the function above the x-axis, when attraction is stronger 
than repulsion, represents the sync phase of our algorithm when attraction is 
the stronger force overall and highly interconnected vertices get relatively close 
to each other in the layout. Effectively, in this phase vertices find their approxi¬ 
mate placement in the final layout. However, it is likely that some vertices will 
be placed too close to each other which makes the outcome of this phase aesthet¬ 
ically unacceptable, thus, we continue the iterations. The part of the function 


below the a;-axis, when the repulsion becomes stronger overall (and grows in¬ 
creasingly faster than the overall attraction), represents the burst phase of our 
algorithm when vertices are spaced out for achieving an aesthetically pleasing 
final result. 



different choices of M{t). 


3.2 Uniform Force Magnitude 

The turning point between the sync and the burst phase is when the total mag¬ 
nitude function /(t) becomes equal to zero, i.e. when the first term 
in Equation ([^ becomes equal to the second term M(t)n{n— 1). If this happens 
for t = tp then the value of M(tp) is the one shown in Equation Q. For a graph 
with n > 2 and m > n — 1 (which is a reasonable assumption) M(tp) 3> 1. 
Thus, by choosing M(l) < M(tp) and having M{t) increase monotonically, the 
algorithm starts with syncing the placement of highly interconnected vertices in 
the first iterations and then, as M{t) becomes larger than M{tp), vertices are 
spaced out in the burst phase. 


Our preliminary experiments with a few graphs, while tuning the parameters 
of the algorithm, have shown evidence that Sync-and-Burst produces aesthet¬ 
ically pleasing results after a relatively short sync phase followed by a longer 
burst phase. Thus, if the number of iterations is fixed in advance then we want 
M(t) to be a function that grows quickly in order to spend fewer iterations in 
the sync phase than in the burst phase. In order to have our algorithm compa¬ 
rable to other classical force-directed algorithms, we want the total number of 
iterations to be 0{n). Clearly, by varying the rate at which M{t) grows we can 








affect how many iterations does the algorithm perform before M (t) reaches the 
value M{tp). 

Let the total number of desirable iterations in the sync phase be sn and the 
total number of desirable iterations in the burst phase be bn where 0 < s < b. 
We want M{t) to be a monotonically increasing function that grows fast. Thus, 
a suitable fit for M(t) would be either a power or an exponential function. 
The choice we made in this work is to have M(t) = (t M(tp) / snY'^ which is 
expressed in terms of n and m in Equation (§. For evaluating Sync-and-Burst 
we tuned the total number of iterations (s + b)n and the value of the parameter 
s by experiment as explained in section 

, / \ ( Itm'^ \ , . 

-TV (5) 

\ sn^[n — Ij / 


3.3 Vertex Coordinates 


Similar to other force-directed algorithms, Sync-and-Burst starts with a random 
initial layout. At each iteration the coordinates of each vertex Vi are updated 
based on the total force of attraction fa{vi,t) and the total force of repulsion 
fr{vi,t) at Vi- Let Oij[t) denote the angle between the x-axis and the straight 
line connecting vertices Vi and Vj at iteration t. Then we compute the x- and 
y-components of fa{vi,t) as shown in Equations (6a) and (6b), respectively. The 
value Qij in Equations (6a) and (6b) is the fj-element of the adjacency matrix 


of the graph. Similarly, the x- and y-components of frivi,t) are computed as 
shown in Equations (7a) and (7b). 


fa(vi,t).x = mM {ajj cos % (t - 1)) 

(6a) 

3¥=i 


faiv„t).y = (oy sin%(t- 1)) 

(6b) 



fr{Vi,t).X = -M(t)^COS%(< - 1) 

(7a) 



fr{vi,t).y = -sin 0ij{t - 1) 

(7b) 



Vi{t).X = faiVi,t).X -f friv„t).X 

(8a) 

+ 

II 

(8b) 


Unlike a typical force-directed algorithm, at each given iteration Sync-and- 
Burst does not move vertex Vi to a new position starting from its previous 
position. Instead, the new position of vertex Vi depends on its previous position 










indirectly. It is influenced by the positions of all vertices relative to Vi at the 
previous iteration through the angles % (see Equations (8a) and (8b)). Finally, 


in order to draw the graph, we normalise the computed coordinates by scaling 
them to the range [0,1]. 


3.4 The Complete Algorithm 

The complete Sync-and-Burst algorithm is presented as Algorithm The body 
of the outermost for- loop in lines [5|20| is a single iteration of Sync-and-Burst. In 
total, we have {s-\-b)n iterations, where s and b are integer constants which repre¬ 
sent the lengths of the sync and burst phases, respectively. The algorithm starts 
with an initial random layout. Since the first iteration of the algorithm needs the 
previous value of the uniform force magnitude, i.e. M(0), we set M(0) = 1/m 
with the intention that we start with a relatively small magnitude below 1. Then 
for t > 1 the value oi M{t) is computed according to the formula in Equation ([^. 

For each vertex Vi, the x- and y-components of the total attraction force 
applied to it are computed in lines [TOp^ and the influence of the total repulsion 
force applied to it is taken into account in lines I4p7 The values dijiO) in the 
first iteration of the outermost for- loop are computed from the initial random 
layout. Finally, we have the coordinates of vertex vi at iteration t in lines [T9]|^ 
If the layout at iteration t has to be drawn, we scale all coordinates to fit in the 
interval [0,1]. 


Algorithm 1 Sync-and-Burst for graph G = {V,E) with \V\ = n and \E\ = m 
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M(0) 
for t < 
M{t) - 
for i i 

fx^ 


place all vertices at randomly picked positions in the drawing area 
s— 1/m 

- 1 to (s -I- b)n do 

{2tm^ / [sv? {n — 1)))^° 

1 to n do 
- 0 
fVi ^ 0 

for / 1 to n do 

if {vi,Vj} G E then 

fxi ■<— fxi + mM(t — 1)°'® cos 6ij(t — 1) 
fVi fVi +mM{t- 1)°-^ sm9ij{t - 1) 

end if 
if i 7 ^ j then 

fxi ■(— fxi — M{t — 1) cos 9ij {t — 1 ) 
fVi fVi - M{t - 1) sin 9ij{t - 1) 


end if 
end for 

Vi{t).X <— 
Vi{t).y -s- 

end for 
end for 


fXi 

fVi 















The computational complexity of Sync-and-Burst is 0{n^) since the out¬ 
ermost loop is executed (s -I- b)n times and its body contains two loops, each 
visiting all vertices. Thus, Sync-and-Burst is comparable to other classical force- 
directed algorithms, including Fruchterman-Reingold, which have the same com¬ 
putational complexity [S]. As reported in section]^ we have experimentally found 
that Sync-and-Burst gives aesthetically pleasing results in 20n iterations. 


4 Experimental Results and Discussion 


We ran both Sync-and-Burst and Fruchterman-Reingold for the Rome Graphs 
introduced by Di Battista et al. [3]. We used all 11534 Rome Graph^which have 
vertex count ranging from 10 to 100 vertices. In order to tune the parameters of 
Sync-and-Burst we ran a pilot study with a few selected graphs. We present the 


layouts of some of them in section 4.2 We discovered that Sync-and-Burst (with 
the aforementioned choice of the M (t) function) achieves aesthetically pleasing 
results in 20n iterations, where n is the number of vertices, and the length of the 
sync phase is inversely proportional to the variance of betweenness centralities in 
the graph. That is, the higher the variance is, the shorter the sync phase needs to 
be. In particular, we achieved aesthetically pleasing results for all graphs in the 
pilot study by setting the parameter s to 20 divided by the standard deviation 
of betweenness centralities in the graph, but keeping s at most 4. That is, the 
sync phase takes at most 20% of all 20n iterations and is typically shorter. 

We created a Flickr picture gallerj0which contains both Sync-and-Burst and 
Fruchterman-Reingold layouts of selected graphs (including some of the Rome 
Graphs). 


4.1 Rome Graphs 


For each layout of a Rome graph we recorded a number of characteristics sum¬ 
marised in Fig.[^ These are all standard graph drawing aesthetic criteria such as 
the number of edge crossings, the average angles between adjacent and crossing 
edges, the standard deviation of edge lengths and a numerical expression of the 
distribution of vertices throughout the drawing area. One way of quantifying the 
distribution of vertices is by finding the minimum distance between two vertices 
in the layout. The higher this value is the more evenly distributed the vertices 
are [9]. We also used a second method of quantifying the distribution of vertices. 
Let d* be the distance between vertex Vi and the vertex which is closest to it in 
the layout; let d** be the distance from Vi to the closest point on the border of 
the drawing area; and let = min{fi*/2, d**}. Then, what we call vertex distri¬ 
bution D in Fig. 2(c) is the ratio of the total area covered by the circles with 
radii and the actual drawing area A as expressed in Equation ([^. The closer 
D is to 1 the more evenly vertices are distributed throughout the drawing area. 


^ available at http://graphdrawmg.org/data.html 
^ https://www.flickr.eom/photos/133882099@N02 





In all our experiments, the drawing area is the smallest rectangle that the lay¬ 
out (either Sync-and-Burst or Fruchterman-Reingold) can fit in, scaled (together 
with the layout) to a rectangle the larger side of which has a unit length. 


D = 


E n 

i=in 


(9) 


For summarising the results, we partitioned the graph dataset into “buckets” 
according to their vertex count, placing a graph with n vertices into bucket [n/5j. 
The values in Fig. are the average for a bucket. The results in Figs. |2(a)] and 
|2(b)| show that there is no clear winner between Sync-and-Burst and Fruchter¬ 
man-Reingold when it comes to the number of edge crossings and angles be¬ 
tween crossing and adjacent edges. While the Fruchterman-Reingold layouts have 
slightly lower number of edge crossings, the Sync-and-Burst layouts tend to have 
slightly larger angles between both crossing and adjacent edges for Rome Graphs 
with vertex count above 55 but the difference is negligible. 

The actual difference between the two algorithms can be clearly seen in 
Figs. 2(c) 2(d)| and |2(e)| Sync-and-Burst clearly distributes vertices increas¬ 
ingly more evenly throughout the drawing area as the size of the graph grows 
(see Figs. 2(c)| and |2(d)[ ). The price that Sync-and-Burst pays for the even distri¬ 
bution of vertices is the larger variance of edge lengths (see Fig. 2(e)). Classical 
force-directed algorithms achieve compact drawings with relatively uniform edge 
length, however, frequently also with large unused (blank) space and entangled 
layouts of dense subgraphs. Sync-and-Burst, on the other hand, achieves much 
better utilisation of the drawing area with more even distribution of the vertices 
by allowing longer edges. The results in Figs. 2(c) and |2(d)| also suggest that 
Sync-and-Burst may scale up better than Fruchterman-Reingold as it tends to 
give increasingly better results in terms of vertex distribution as the size of the 
graph grows. 

Finally, we have compared the running times of Fruchterman-Reingold and 
Sync-and-Burst on a desktop computer with 8 GB RAM and a 3.20 GHz quad- 
core CPU. Fig. |2(f)] shows that an iteration of Sync-and-Burst is faster than an 
iteration of Fruchterman-Reingold. This is something we expected because Sync- 
and-Burst uses simpler force magnitudes which do not depend on the distance 
between vertices and thus are faster to compute. 


4.2 Selected Graphs 

In this section we demonstrate the work of Sync-and-Burst on a few selected 
graphs which we also used in a pilot experimental study to tune parameters of the 
algorithms. The layouts of 6 of these graphs are shown Figs.[3][^ For comparison, 
we also present a Fruchterman-Reingold layout of each of these graphs. In the 
caption of each layout we denote Sync-and-Burst by SnB and Fruchterman- 
Reingold by FR. The vertex distribution, the number of edge crossings and the 
standard deviation of edge lengths for the layouts in Figs. |3][^ are presented in 
Table [l] We counted an edge crossing for each pair of edges which cross, thus in 































(e) 


(f) 


Fig. 2. Comparison between Sync-and-Burst and Fruchterman-Reingold layouts of the 
Rome Graphs: number of edge crossings (a), average angle between any pair of adjacent 
edges and average crossing angle; the average crossing angle for planar layouts was 
set to 90° (b), vertex distribution (c), minimum distance between a pair of vertices 
multiplied by the number of vertices (d), standard deviation of edge lengths (d), and 
running time in seconds (d). 












the Sync-and-Burst layout of the Wagner graph (see Fig. |5(a)| ), for example, we 
counted 6 edge crossings instead of 1. 

The layouts of the selected graphs confirm the findings in section 4.1 that 
Sync-and-Burst distributes vertices evenly throughout the drawing area. Some¬ 
times, this is not in terms of the vertex distribution we have measured according 
to Equation but instead in terms of producing a highly symmetrical layout. 
Such are the symmetrical Sync-and-Burst layouts of the 15 x 5 queen graph 
(see Fig. |3(c)|, the Wagner graph (see Fig. |5(a)|) and the Heawood graph (see 
Fig.[5^ 3 whose Fruchterman-Reingold layouts have higher (i.e. better) vertex 
distribution according to Equation 0 (see Table [^. Our conjecture is that 
by producing circular-shape layouts and allowing long edges, Sync-and-Burst 
is able to highlight symmetries in graphs which do have symmetrical layouts 
when relatively long edges are permitted. The layouts of the selected graphs also 
demonstrate that Sync-and-Burst achieves good results for dense and scale-free 
graphs (see Figs. 3(a)[ 3(c)| and 4(c) I which further suggests that it may scale 
up better than Fruchterman-Reingold. 






(a) 8 X 8 SnB 


(b) 8 X 8 FR 


(c) 15 X 5 SnB 


(d) 15 X 5 FR 


Fig. 3. Sync-and-Burst layouts (a) and (c) compared to Fruchterman-Reingold layouts 
(b) and (d) of the 8x8 and the 15 x 5 queen graphs. Sync-and-Burst achieves layouts 
which are clearly more symmetrical. 





Fig. 4. Sync-and-Burst layouts (a) and (c) compared to Fruchterman-Reingold layouts 
(b) and (d) of a sparse disconnected graph (graphl) and a randomly generated scale- 
free graph with 130 vertices (graph2). 






















(a) Wagner SnB (b) Wagner FR (c) Heawood SnB (d) Heawood FR 


Fig. 5. Sync-and-Burst layouts (a) and (c) compared to Fruchterman-Reingold layouts 
(b) and (d) of a the Wagner graph and the Heawood graph, respectively [TT] . 


Table 1. Vertex distribution, number of edge crossings and standard deviation of edge 
lengths for the layouts in Figs. |3|5| 


graph 

n 

m 

vertex 

distribution 

number of 
edge crossings 

stdev of 
edge lengths 

SnB 

FR 

SnB 

FR 

SnB 

FR 

queen 8x8 

64 

728 

0.480 

0.409 

31322 

30462 

0.469 

0.358 

queen 15 x 5 

75 

935 

0.088 

0.391 

41950 

44064 

0.521 

0.297 

graphl 

24 

32 

0.142 

0.004 

1 

1 

0.095 

0.009 

graph2 

130 

190 

0.299 

0.240 

287 

308 

0.281 

0.070 

Wagner 

8 

12 

0.066 

0.118 

6 (or 1) 

5 

0.582 

0.226 

Heawood 

14 

21 

0.038 

0.217 

14 

13 

0.605 

0.149 


5 Conclusions 

Sync-and-Burst, the algorithm presented in this paper, falls into the category 
of classical force-directed algorithms. While following the general schema of the 
algorithm of Fruchterman and Reingold, Sync-and-Burst uses forces of attrac¬ 
tion and repulsion whose magnitude is independent from the distance between 
vertices. Instead, a uniform magnitude of attraction and a uniform magnitude 
of repulsion are applied throughout the graph at each iteration and these mag¬ 
nitudes monotonically increase as the number of iterations grows. As a result of 
this, the Sync-and-Burst layouts are always circular in shape. This could be an 
advantage for applications where users prefer to see a variety of graphs drawn 
in the same fashion. Moreover, the experimental study presented in section 
demonstrates that Sync-and-Burst distributes vertices more evenly throughout 
the drawing area than the algorithm of Fruchterman and Reingold, while hav¬ 
ing similar number of edge crossings and angles between crossing and adjacent 
edges. Thus, we conclude that Sync-and-Burst distributes vertices more evenly 












while still achieving a clear and aesthetically pleasing layout. In particular, we 
achieved very good results for graphs with symmetries in their structure. 

Future work on this algorithm may potentially reveal a better choice for the 
magnitude function M(t) and the number of iterations, as well as for the ratio 
between the lengths of the sync and burst phases. In either case, because of 
the simplicity of the forces, an iteration of Sync-and-Burst can be faster than an 
iteration of other classical force-directed algorithm (as reported in section]^ and 
we have shown that Sync-and-Burst layouts are aesthetically pleasing in 0{n) 
iterations. 
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